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(同济 大 学 环境 科学 与 工程 学 院 ， 上 海 200092) 
摘 要 垃圾 焚烧 会 产生 大 量 的 二 次 污染 , 使 得 原本 极为 推崇 的 垃圾 焚烧 发 电 倍 受 质疑 。 本文 利 用 CFD 技术 对 垃圾 焚烧 产生 的 
烟 气 通过 烟 轴 高 空 扩 散 的 规律 进行 数值 模拟 ， 研 究 在 多 尺度 空间 条 件 下 流 场 和 污染 物 分 布 。 采 用 标准 -e 消 流 模型 及 组 分 输 
运 模型 ， 分 析 不 同 排 烟 高 度 和 速度 对 排 烟 地 点 下 风向 1.4 公里 小 区 建筑 群 内 污染 物 分 布 的 影响 。 研 究 表明 SO，、 二 恶 英 等 烟 气 
污染 物 在 小 区 建筑 间 的 浓度 分 布 随 着 排 烟 高 度 和 速度 的 增加 而 减弱 , 但 由 于 气流 分 布 产 生 涡 旋 的 影响 使 得 建筑 物 背 风 面 污染 物 
浓度 较 高 。 本 文 将 计算 结果 与 正 态 分 布 假设 下 的 高 斯 烟 羽 模型 计算 结果 进行 比较 ， 两 者 有 较 好 的 一 致 性 。 
关键 词 垃圾 焚烧 ， 数 值 分 析 ; 二 恶 英 ;高 斯 烟 羽 模型 
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Multiscale Simulation and Analysis on the Diffusion of 


the Waste Incineration Flue Gas 


CUI Peng-Yi WU Zhi-Gen Tao Wen-Quan* 

(College of Environmental Science and Engineering, Tongji University, Shanghai 200092, China) 
Abstract The large number of secondary pollutants, which were produced from waste incineration flue gas, would 
cause the highly admired waste incineration power generating project into questions. The CFD technology was 
applied to simulate the diffusion characteristics of waste incineration flue gas from chimney, to investigate the flow 
and pollutant distribution at multiscale space conditions. The standard k - ¢ turbulence model and component 
transport model were adopted to evaluate the pollutant distribution influence between buildings at downwind 1.4km 
from discharge position with different exhaust height and velocity. It was found that the concentration distribution of 
flue gas pollutant, such as SO2 and TCDD, decreased with the increase of exhaust height and velocity, and a higher 
pollutant concentration on the leeward of buildings was caused by the vortex induced from wind field. Finally, the 
analysis results of Gaussian plume model were found to be in a good agreement with the simulation results. 
Key words waste incineration; numerical simulation; TCDD; Gaussian plume model 


0 引言 

城市 生活 垃圾 的 主要 处 理 方式 有 卫生 填 埋 、 堆 最 著名 的 是 高 斯 烟 羽 扩散 模型 00, 实验 方法 是 一 类 
肥 和 焚烧 。 其 中 ， 近 30 年 中 来 发 展 起 来 的 垃圾 焚 依据 观测 资料 的 分 析 方法 ,其 精度 很 大 程度 上 取决 
烧 技术 是 目前 较为 经 济 有 效 的 垃圾 处 理 技术 之 一 于 观测 密度 。CFD 技术 是 近年 来 迅速 发 展 起 来 的 一 
03 习 。 然 而 ， 垃 圾 焚烧 会 产生 大 量 的 二 次 污染 ， 所 种 研究 流动 传 质 、 传 热 与 反应 等 的 有 效 手段 ， 在 环 
生成 的 烟 气 中 含有 至 今 毒性 最 强 的 物质 二 恶 英 外 。 境 工程 领域 作为 一 种 研究 手段 越 来 越 受 到 人 们 重 
目前 , 全 球 二 恶 英 的 产生 量 几乎 全 部 来 自生 活 垃圾 视 02431。 与 实验 方法 相 比 , CED 技术 具有 经 济 、 高 
焚烧 厂 排放 出 来 烟 气 回 。 效 的 特点 ， 相 比 于 理论 分 析 ， 其 能 够 求解 更 为 复杂 

烟 气 从 烟 身 排出 后 扩散 规律 的 研究 主要 包括 的 流动 扩散 问题 。 本 研究 利用 FLUENT 软件 ， 构 
理论 分 析 、 实 验方 法 和 数值 模拟 个 中 。 理 论 分 析 中 建 多 尺度 数值 计算 模型 ,研究 垃圾 焚烧 产生 的 烟 气 
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经 不 同 高 度 和 速度 排放 之 后 , 在 大 气 中 扩散 传 质 的 
规律 , 对 大 气 环境 污染 控制 与 评价 具有 一 定 指导 意 
义 。 


1 数学 模型 及 求解 


1.1 几何 模型 与 网 格 划 分 

如 图 1 所 示 , 计算 区 域 为 800 m x 800 m x 3000 
m 长 方 体 区 域 , 前 后 两 个 正方 形 分 别 是 计算 区 域 的 
ADAH. WAIN EE 5m WEE, A A, E 
计算 中 五 取 值 为 25 m, 50 m 和 100 m; 烟 向 位 于 
入 口 平 面 下 风向 400 m Xb, 小 区 位 于 入 口 平 面 下 游 


1.8 km 处 ， 也 就 是 烟 向 下 游 1.4 km 处 ， 小 区 内 建 

筑 的 尺寸 如 图 1 中 放大 部 分 所 示 。 计算 区 域 的 选择 

考虑 到 烟 奥 排出 的 烟 羽 有 充分 发 展 的 空间 , 不 会 受 
到 计算 空间 四 周边 界 的 影响 , 保证 计算 区 域 的 出 口 
可 以 采用 充分 发 展 流动 的 边界 条 件 。 


图 1 计算 区 域 


Fig.l Computational domain 


本 模型 计算 区 域 较 大 , 流动 的 特征 尺度 从 整个 
计算 区 域 的 几 千 米 到 住宅 小 区 的 几 十 米 , 属于 系统 
多 尺度 问题 (“站 。 系 统 多 尺度 问题 可 以 采用 统一 方 
程 整 场 求解 以 及 分 层 求解 逐 级 细 化 两 种 方法 , 根据 
本 文 的 具体 情况 采用 整 场 求解 的 方法 , 在 网 格 的 划 
分 上 充分 考虑 了 不 同 区 域 对 网 格 朴 密 的 要 求 和 计 
算 精 度 与 求解 速度 之 间 的 平衡 。 网 格 尺度 从 一 米 到 
te 

染 物 影响 的 小 区 周围 以 及 排放 源 到 小 区 之 间 污 
feet 采用 较 高 的 网 格 密度 ; 在 其 他 
污染 物 很 难 到 达 的 高 空 及 烟 羽 影响 较 小 的 两 侧 采 
用 较 粗 的 网 格 。 如 图 2 所 示 ， 烟 秽 附近 采用 0.1-0.5 
米 的 分 辨 紊 的 网 格 尺度 , 图 3 所 示 的 小 区 及 附近 采 
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用 0.5-1 米 的 网 格 尺 度 ， 而 在 其 两 侧 及 500 KUE 
的 高 空 网 格 尺度 则 在 十 到 几 十 米 不 等 。 总 体 网 格 数 
在 400 万 左右 ， 网 格 利 用 Gambit 软件 生成 。 


图 2 AA el Ae BASEL PX 
Fig.2 Chimney and nearby grid 


图 3 小 区 及 附近 网 格 
Fig.3 Community and nearby grid 


1.2 控制 方程 与 边界 条 件 

将 烟 气 排放 的 扩散 传 质 过 程 视 为 不 可 压缩 流 
体 流动 过 程 ， 由 于 烟 气 排放 温度 较 高 ， 对 其 密度 影 
响 较 大 , 一 般 排放 口 的 烟 气 密度 p, 与 环境 大 气 的 密 
FE po 不 同 , 在 垂直 方向 的 动量 方程 中 要 考虑 到 浮 升 
力 效应 的 影响 。 在 其 余 方程 中 ， 可 以 用 环境 大 气 密 
度 po 代 营 烟 气 密度 p。， 控制 方程 由 连续 性 方程 、 动 
量 方程 、 能 量 方程 、 组 分 输 运 方程 以 及 封闭 上 述 方 
各 的 济 流 方 程 即 标准 . e MARAR, 控制 
方程 由 商业 软件 FLUENT 求解 。 

FLUENT 软件 中 参数 的 设 定 : 满 流 施 密 特 数 
Sct 取 值 0.7; 系数 C, HUH 0.0845; 本 文 不 考虑 烟 气 
压力 、 温度 以 及 其 他 组 分 变化 对 SO 和 二 恶 英 扩散 
系数 的 影响 , 将 扩散 系数 视 为 常数 , 取 为 2.88 x105 
msl， 计算 区 域 入 口 采用 速度 入 口 边界 条 件 ， 大 
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气 边 界 层 速度 入 口 边 界 符 合 以 下 方程 : 1.80E.015 

U/Un=ln(z/zo)/In(2w/Z0)s 其 中 , z 为 距离 地 面 的 高 度 ; 60E.016 LE 

20 为 地 面 粗糙 高 度 ， 设 为 0.3 m; za AARP 140E.015 上 o 

的 高 度 ; Un 为 烟 向 排放 口 所 在 高 度 环境 风速 ， 设 3 1.206-016 L 

为 4.5 m/st' 1。 入 口 环境 大 气温 度 为 288 K， 烟 气 等 1007-015 F 

污染 物质 量 浓度 为 0。 烟 气 的 入 口 设置 在 烟 奥 的 底 3 8.00E-016 | 

部 ， 该 处 采用 速度 入 口 条 件 ， 该 入 口 设置 烟 气温 度 l1 6.00E-016 |- 

Fy 433 K, SO, 的 质量 浓度 为 8.99x107, Z WHH 4.00E-016 L 

质量 浓度 为 8.33x10". AKERA OREA N H 2.00E-016 上 | TN 
WH, WHAA Op/On=0, p 代表 所 有 求解 变 0.00E+000 L - = ———- 
Eo ART APA RAO, RR JA i BE /m 

界 条 件 要 求 垂 直 该 边界 的 速度 分 量 的 值 为 0， 而 其 

它 标量 及 和 平行 于 该 边界 的 速度 分 量 的 一 阶 导 数 为 5V = 6 m/s 下 二 恶 英 平均 质量 浓度 
0。 地 面 、 小 区 建筑 和 烟 秽 表面 为 速度 无 滑 移 ， 物 Fig.5 Average mass fraction of TCDD at V=6 m/s 


质 无 渗透 的 壁面 边界 。 
2.2 排 烟 速 度 对 烟 气 扩散 的 影响 


TE 图 6 给 出 了 在 相同 排 烟 高 度 (50 my) 不 同 排 烟 
2 结果 分 讨论 
m= 6 琳 分 析 及 讨 速度 条 件 下 , 小 区 建筑 范围 内 SO 和 二 恶 英 在 人 呼 
TO 本文 在 相同 大气 条 件 下 , 对 垃 航 焚烧 烟 气 在 相 吸 截面 (z= 1.5 m 所 达到 的 最 大 质量 浓度 。 排 烟 
Do 同 排 烟 温度 不 同 排 烟 高 度 和 排 烟 速度 条 件 下 ,研究 速度 的 增加 意味 着 单位 时 间 污染 物 释放 量 增加 , 各 
了 SO, 和 二 亚 英 等 污染 物 在 大 气 中 扩散 的 规律 以 ste Yc: Shun eet ike RE ABE es ins 
及 对 下 风向 1.4 km 处 居民 小 区 的 影响 ,并 做 了 细致 点 污染 物 的 浓度 值 也 应 随 之 成 比例 增加 ; 但 排 烟 速 
e i Í 度 的 增加 会 提高 烟 气 上 升 的 动量 ， 烟 气 的 抬升 高 度 
2.1 排 烟 高 度 对 烟 气 扩散 的 影响 AH 和 有 效 源 高 He (He=H+AH) 也 将 随 之 增加 ， 
= 4、 图 5 表示 在 相同 排 烟 速度 (6 m/s) 条 件 根据 高 斯 烟 羽 模型 理论 号 有 效 源 高 的 增加 有 利于 
F, SO, 和 二 恶 英 在 z=1.5 m 和 z=24m 小 区 截面 烟 气 中 污染 物 的 高 空 扩散 ， 从 而 降低 其 对 近 地 面 的 
> i i ent L 污染 程度 。 因此 排 烟 速 度 的 增加 对 近 地 面 污染 程度 
S 增加, 小 区 居民 生活 范围 内 SO 和 一 恶 英 浓度 明显 有 两 面 性 .根据 图 6 本 文 所 建 模型 的 情况 分 析 可 知 ， 
“全 降低 ， 排 烟 高 度 从 25 m、50 m 到 100 m， 两 种 污 烟 气 的 排放 速度 由 6 m/s. 10 m/s 增加 到 15 m/s 
— 染 物 质量 请 度 都 有 近 一 个 数量 级 的 减少 , 可见 排 烟 ee n 
= 高度 的 增加 可 以 在 很 大 程度 上 使 得 污染 稀释 扩散 SO: ARRAS A EAT X RE 
到 高 空 大 气 ， 从 而 大 大 减弱 对 下 游 近 地 面 附近 的 污 烟 气 排放 速度 的 增加 而 减弱 , 说 明了 高 空 大 气 对 烟 
染 程度 。 气 的 稀释 程度 要 大 于 烟 气 排放 量 的 增加 对 近 地 面 
2.00E-008 
1.80E-008 | m 350E009 py a — 3.50E-016 
1.60E-008 L 一 @ 一 截面 z= 1.5 m m 3.00E-009 上 ®, lease. 
re —m— 截面 z= 24 m x = 
ual ne Se J 250E-016 = 
= 1.20E-008 L K Aaea = 
et [ ~ 4 2.00E-016 | | 
ux 1.00E-008 | 种 150E-009 上 ha 
$F 3.00E-009 上 z nell: ~~ 了 iene 
sans = 4 A =» 
4.00E-009 L Re E e = 
.00 ee | a celia Te e a 0 a a a Mn e Y 
F @—— E 7 
0.00E+000 = = a = 排 烟 速度 (m/s) 
DANE 图 6 不 同 排 烟 速度 下 SO 和 二 恶 英 最 大 质量 浓度 
Al4V = 6 m/s 下 SO; 平均 质量 浓度 


Fig.6 Maximum mass fraction of SO, and TCDD at 


Fig.4 Average mass fraction of SO, at V=6 m/s different exhaust velocity 


的 影响 程度 ,也 进一步 说 明 提高 烟 气 排放 的 有 效 源 
高 〈 包 括 增加 烟 向 的 几何 高 度 ) 有 利于 减弱 烟 气 污 
染 物 对 近 地 面 的 影响 。 

2.3 楼 宇 间 的 街区 峡谷 对 烟 气 扩散 的 影响 

图 7 显示 了 楼 宇 间 街区 中 SO, 质量 浓度 分 布 的 
模拟 结果 。 由 图 可 见 ,小 区 建筑 间 与 风 主 方向 平行 
的 街区 SO, 等 污染 物 浓度 要 明显 低 于 与 来 风 主 导 
方向 垂直 即 建筑 背风 面 的 街区 ,这 种 现象 与 建筑 物 
布置 对 环境 风速 度 场 的 影响 有 关 , 从 图 8 可 以 看 出 ， 
建筑 的 背风 街区 产生 明显 的 涡流 , 不 利于 新 风量 的 
更 新 ， 也 不 利于 污染 物 的 随 风 扩 散 稀 释 ， 造 成 污染 
物 浓度 积累 ,使 得 建筑 背风 街区 局 部 空气 质量 明显 
下 降 。 
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图 7 SO, 质量 浓度 在 z=1.5m 截面 分 布 


Fig.7 Mass fraction on the plan at z=1.5m 
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图 8 街区 竖 直 截面 速度 分 


Fig.8 Velocity vector distribution on vertical plan 


2.4 模拟 结果 与 烟 羽 模型 理论 分 析 结果 比较 

在 大 气 环 境 流 场 稳 定 ， 忽 略 横向 、 竖 向 流速 及 
纵向 弥散 作用 ， 对 于 有 效 排放 高 度 为 He 的 连续 排 
放 点 源 , 污染 物 浓度 的 空间 分 布 符合 连续 排放 点 源 
的 高 斯 模式 中 。 如 图 9 所 示 ， 本 模型 轴线 y = 400 
m 近 地 面 $S0; 质 量 浓度 分 布 与 高 斯 烟 羽 模型 理论 值 
有 较 好 一 臻 性。 本文 在 污染 排放 几何 高 度 为 50 m, 
排放 速度 10 m/s, WHEE 433 K 等 条 件 下 ， 地 面 (z = 
0 m)SO> 质量 浓度 在 2 km 左右 达到 最 大 , 理论 分 析 
在 1.8 km 左右 为 最 大 。 由 于 地 面 轴线 在 x= 1.8 km 


处 受 小 区 建筑 群 的 影响 , 使 得 本 模型 近 地 面 污染 物 
浓度 模拟 与 理论 计算 值 之 间 存 在 一 定 差别 。 


B 
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图 9 SO, 浓度 分 布 理 论 与 模拟 值 比 较 
Fig.9 Comparison between theoretical and simulation 
values of SO, concentration distribution 


3 讨论 

本 文中 应 用 计算 流体 力学 技术 对 生活 垃圾 焚 
烧 产 生 的 烟 气 污染 物 通 过 烟 秽 高 空 排放 后 ,在 大 气 
环境 中 的 稀释 扩散 过 程 和 规律 进行 了 多 尺度 数值 
模拟 研究 , 计算 在 包含 污染 源 排放 点 及 其 下 风向 影 
响 小 区 在 内 的 一 个 足够 大 的 空间 进行 , TS TAA 
高 度 及 排 烟 速度 的 影响 。 分 析 了 不 同 条 件 下 污染 源 
(SO;、 二 恶 英 等 污染 物 ) 在 空间 扩散 的 规律 及 在 
小 区 内 分 布 特点 , 对 于 垃圾 焚烧 厂 附近 环境 影响 评 
价 有 一 定 的 指导 意义 。 本 模型 考虑 的 是 在 稳定 的 气 
象 条 件 下 ， 污 染 物 随 气流 的 扩散 ， 但 大 气 环境 瞬息 
万 变 ， 局 部 气流 变化 具有 随机 性 ， 因 而 对 大 气 中 污 
染 物 的 扩散 有 很 大 的 影响 ,也 就 是 说 烟 气 等 大 气 污 
染 物 排放 到 大 气 中 后 , 其 运动 扩散 规律 有 很 大 随机 
性 ， 如 何 把 握 这 种 规律 也 是 以 后 研究 的 重点 。 


4 结论 

本 文 利用 CFD 技术 在 稳 态 条 件 下 对 垃圾 焚烧 
产生 的 烟 气 污染 通过 烟 肉 扩散 的 规律 进行 了 多 尺 
ERUR, KA MIRE k- e 两 方程 模型 及 组 分 
输 运 模型 ， 研 究 了 在 较 大 区 域 范围 条 件 下 ， 烟 气 在 
闻 上 气流 分 布 和 污染 物 分 布 ; 分 析 了 不 同 排 烟 高 
和 排 烟 速度 条 件 下 对 距离 排 烟 地 点 1.4 公里 处 的 
烟 气 污染 分 布 的 影响 。 研 究 表明 SO,、 二 恶 英 
气 污染 物 在 小 区 中 的 浓度 分 布 随 着 排 烟 高 度 
烟 速 度 的 增加 而 减弱 ， 由 于 气流 分 布 产生 涡 旋 
J 影响 使 得 建筑 物 之 间 的 背风 街区 污染 物 浓度 较 
高 。 计算 结果 与 正 态 分 布 假设 下 的 高 斯 烟 羽 模型 进 
行 比 较 ， 两 者 有 较 好 的 一 致 性 。 
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